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We present the results for hadronic spectrum, decay constants and the light quark masses obtained with 
non-perturbatively improved Wilson fermions. We also give our preliminary results for the heavy-light decay 
constants. 
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The method for systematic improvement of the 
Wilson action and bilinear quark operators, re- 
cently proposed by ALPHA-collaboration 0] , pro- 
vides a full cancellation of the linear discretization 
errors. In other words, the results of lattice cal- 
culations of any on — shell quantity, do converge 
to the continuum limit (a — > 0) much faster , i.e. 
with the rate 0(a 2 ). 

To verify and quantify the effect of this im- 
provement (first improved studies of the spec- 
troscopy were done in ||), we generated 100 inde- 
pendent gauge field configurations on a 24 3 x 64 
lattice, at (3 = 6.2, in the quenched approxima- 
tion. We used the non-perturbative(NP) value 
of the coefficient for the improvement of the ac- 
tion, c sw , which guarantees the full improvement 
of the hadron mass spectrum. The values of all 
other counterterms coefficients, as well as of the 
renormalization constants which we used in this 
work, are listed in Tab. [I]. Note that the NP 
results for Zs, Zp and Zt (given in Tab. [j] at 
[ia, = 1), are new. We calculated several correla- 
tors for mesons made of four degenerate quark 
'flavors'. In the analysis of all hadronic quanti- 
ties (spectrum and decay constants), we used the 
method of 'physical lattice planes', proposed and 
extensively used in Ref. B. Further details on 
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Table 1 

Renormalization constants Zj, and the improvement co- 
efficients cj and bj at /3 = 6.2 (c sw = 1.614). For the 
constants which are not computed non-perturbatively (x), 
we take the result of the 1-loop boosted perturbation the- 
ory (BPT). 
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the computation and the analysis can be found 
in Here, we will address several issues which 
are relevant for the improvement, and report our 
results. 

Critical hopping parameter: The first observa- 
tion that we would like to make, is that the value 
K cr n must be reached through a quadratic fit: 

ax(2m q ) + a 2 (2m q ) 2 (1) 



M PS 



where 2m q = k — n crit , is the quark mass 



that can be derived from the vector Ward iden- 
tity (VWI). Namely, in order to be consistent 
with the O(a)-improvement, the usual (linear) 
fit in quark mass should also be modified as 
Mp s ~ rhq, where rh q = m q (l + b m m q ). In 
perturbation theory, the coefficient b m < 0. By 
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fitting; our data with (@), we unambiguously ob- 
tain a-i > 0, which means that the positive cur- 
vature, due to the physical m q corrections, ex- 
ceeds the lattice artifacts correction of the same 
order (m^). The values of K cr n obtained from 
the quadratic and the linear fits differ appreciably 
( K q c Z dllm] = 0.1 3585(2) [O.f 3576(2)]). One more 
convincing argument in favor of K^if, is provided 
by the O(a) improved axial Ward identity (AWI), 
by which the quark mass is defined as: 



Pi 



2<P(i)0t(O)) 



G(a 2 ). 



(2) 



From the linear fit p q (l/n), we get lim P(! ^o K = 
0.13584(5). A quadratic fit p q (l/n) gives the 
same result and we deduce that t\. r ^ t is the cor- 
rect value. 

Calibration of a -1 : We adopt the argument (of 
Ref. H) that, due to the minimal extrapolation 
required, the most convenient quantity to fix a -1 
is m K ,. We obtain: a -1 = 2.75(17) GeV . An en- 
couraging feature of the improvement is that the 
value of the inverse lattice spacing, as obtained 
by using various physical quantities (e.g. m p , fx, 
/tt), remains stable and close to the one obtained 
from the computation of the string tension (the 
pure gluonic quantity). 

Light quark masses: The quark mass can be ex- 
tracted in two ways, depending on the defini- 
tion we chose, VWI or AWI. The two bare quark 
masses, when appropriately renormalized, should 
converge to the same physical result. Recently 
in Ref. ||, this was proved to be the case when 
the renormalization constants (RC) are computed 
non-perturbatively. We have also calculated the 
RC's (in the RI scheme) by using the method 
for NP renormalization || (see Tab. 0). Sub- 
tleties and computational details on the extrac- 
tion of bare lattice quark masses can be found 
in 0. The ^-coefficients, which are relevant for 
the improvement out of the chiral limit only, i.e. 
Zs,P.A{m,[i) ~ (I + b s ,p,Ama)Z s ,p,A(0,p), are 
taken from the (boosted) perturbation theory. 
Our final results read: 



4-5(4) w2r , n [4.9(4), 



mT S \2GeV) = 111(12) _ [121(13), 



jMeV, 



Note that the matching between the MS and 
RI schemes was also performed at NNLO accu- 
racy ||, which is the new result. To the same ac- 
curacy, the renorm. group invariant quark masses 
are: 

m uA = 7.2(6) MeV, m s = 177(19) MeV. 

Light meson decay constants: For the full O(a) 
improvement of the quark bilinear operators, the 
coefficients cj, of the appropriate counterterms, 
are tuned to their non-perturbative values. Only 
Ct is to be taken from the boosted perturbation 
theory. By using the improved currents matrix 
elements, we define decay constants in the follow- 
ing way: 
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where V, V denote the pseudoscalar and vector 
mesons at rest. The corrections aFj \ correspond 
to the matrix element of: the gradient of the pseu- 
doscalar density, the divergence of tensor current, 
and to the curl of the vector current, respectively. 
Empirically, aF^p is large, whereas aF^p and 

T(l) 

aF v are small. However, when multiplied by 
corresponding cj, their respective contributions 
never exceed ~ 5%, relative to the leading term 

( ^ 0)) - 

All the bare lattice local currents in (H), 
should be multiplied by appropriate renormaliza- 
tion constants. In the chiral limit, they are all 
computed nonperturbatively (see Tab. [I]). This 
is not the case for bj's, which are all (except for 
by) taken from 1-loop boosted perturbation the- 
ory. With all these ingredients, we calculated the 
pseudoscalar and vector decay constants. The 
results, obtained using a _1 (m^*), are listed in 
Tab. 2. 

It is well known that the ratio fx / f-n — 1 (which 
is the 'measure' of the SU(3) flavor breaking), as 
extracted from the lattice calculation, is always 
smaller than the experimental value (0.22). After 
a careful analysis of all the data (unimproved and 
improved), the lesson we draw is: (fx — fn)/ fn 
from the lattice quenched studies is insensitive 
to the improvement and remains more than 30% 
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Table 2 

Decay constants. 
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lower than its experimental value. 
On the other hand, when compared to the unim- 
proved results, the values for the vector decay 
constants are in substantially better agreement 
with the experimental ones. For comparison 
with the unimproved results, we show the di- 
rectly extracted data (with no extrapolations), 
as obtained from simulations with unimproved, 
tree-level improved and NP improved Wilson 
fcrmions. The other source of vector mesons is 
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Figure 1 . Comparison of directly extracted vector meson 
decay constant with unimproved, tree-level improved and 
NP improved Wilson fermions, at the same j3 = 6.2, and 
on the same lattice, 24 2 X 64. 

provided by the tensor current. The masses ex- 
tracted from vector and tensor correlation func- 
tions are (of course) the same, while the couplings 
are different. In Tab. 2, we give the first lattice 
predictions for fJ K , l p(2GeV). 

Heavy-light meson decay constant (preliminary 
results): In the simulations, we included four 
heavy quarks which we combined with the 
four light ones. After a quadratic extrapola- 
tion/interpolation in light quark masses, we in- 



terpolate in the charm region to estimate fjj and 
fo s • This results in: 

f D = 202 ± 14t° 2 MeV, 

f D , = 231 ±11^ MeV, 

Sd.ISd = 1.11 ±0.03. (4) 

The stability of this result (see M), does not hold 
for B mesons. As usual, we have to extrapolate 
our results by using the heavy quark scaling law: 



fQq ~ 
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(5) 



V m Q? \ mQq 
It turns out that fB,B s are sensitive to whether 
we include the ~ l/ m Q g term in the fit or not. 
This is illustrated in Fig. 2, and its effects, along 
with the scale fixing dependence, are included in 
the systematic error. Our results are: 
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Figure 2. Directly extracted heavy-light pseudoscalar 
decay constants vs. two forms of fit needed to extrapolate 
to the B meson sector. 
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